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Abstract 



o 



We suggest that the combined effect of screening, gluon-induced dissociation, coUisional damping, and reduced feed-down explains 
most of the sequential suppression of T(nS) states that has been observed in PbPb relative to pp collisions at y/sJ/N - 2.76 TeV. 
The suppression is thus a clear, albeit indirect, indication for the presence of a qgp. 
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The suppression of quarkonium states is a very promising 
probe for the properties of the quark-gluon plasma (QGP) that 
is generated in heavy-ion collisions at high relativistic energies. 
In the QGP the confining potential of heavy quarkonium states 
is screened due to the interaction of the heavy quark and the 
antiquark with medium partons and hence, charmonium and 
bottomium states successively melt 1 1 1 at sufficiently high tem- 
peratures Tdiss beyond the critical value Tc - 170 MeV. How- 
ever, additional processes such as gluon-induced dissociation, 
and collisional damping contribute to the suppression. Here we 
concentrate on such processes. 

Charmonium suppression has been studied since 1986 in 
great detail both theoretically, and experimentally at energies 
reached at SPS and RHIC EEglia, and at the LHC Ed. The 
precise origin is still under investigation. Bottomium suppres- 
sion is expected to be a cleaner probe. The T(15) ground state 
with invariant mass 9.46 GeV is strongly bound, the threshold 
to BB decay is at 1.098 GeV. Its lifetime of 1.22-10"2"s is about 
1.7 times as large as the one of 7/^(15) in elementary colli- 
sions. It melts as the last quarkonium in the QGP (depending on 
the potential) only at 4.10 |8 |, whereas the 25 (10.02 GeV) 
and 3S (10.36 GeV) states melt at about 1.6 and 1.2 T^, respec- 
tively. Even at LHC energies the number of bottom quarks in 
the QGP remains small such that statistical regeneration of the 
T states is unimportant. 

T suppression in heavy-ion collisions has recently been ob- 
served for the first time both by the STAR experiment at RHIC 
121, and in 2010 by the Compact Muon Solenoid (CMS) experi- 
ment at LHC ifTOlfTTIl . The latter includes an observation of the 
enhanced suppression of the 2S + 3S relative to the IS ground 
state. CMS data from the 201 1 run shown at this conference 
lfT2ll have much better statistics such that the 2S state can be 
resolved individually. 

In this work we investigate the suppression of 
Y(15), (25), (35) states at LHC energies due to screen- 
ing, gluon-induced dissociation ||T31 . collisional damping, 
and reduced feed-down from the T(25,35) and Xbi^P) and 



XbC^P) states. Whereas gluodissociation below T^. is not 
possible due to confinement, it does occur above where the 
color-octet state of a free quark and antiquark can propagate 
in the medium. The process is relevant below the dissociation 
temperature T^iss that is due to Debye screening, and its signif- 
icance increases substantially with the rising gluon density at 
LHC energies. 

In the midrapidity range \y\ < 2.4 where the CMS measure- 
ment [11 1 has been performed, the temperature and hence, the 
thermal gluon density is high, and causes a rapid dissociation 
in particular of the 25 and 35 states, but also of the 15 ground 
state. At larger rapidities up to the beam value of yheam - 7.99 
and correspondingly small scattering angles where the valence- 
quark density is high lfT4l . nonthermal processes would be more 
important than in the midrapidity region that we are investigat- 
ing here. Thermal gluons will also dissociate the;^'fc(lP) and 
Xbi^P) states which partially feed the T(15) ground state in el- 
ementary collisions fTS'l. 

Due to the small velocity v <sc c of the quarks in the bound 
state, the proper equation of motion for single-particle quarko- 
nium states is the Schrodinger equation, with the color-singlet 
QQ quarkonium potential Vqq. Reasonable parametrizations of 
the potential exist that have been tested in detailed calculations 
of the excited states. 

In particular, the Cornell potential |[T6l has string and 
Coulomb part Vqq = crr-al/r, where cr 0.192 GeV^ (TTI is 
the string tension, and ^ 0.37 the strong- coupling constant 
at the soft scale m^a, (nih - 4.77 GeV) that accounts for the 
short-range gluon exchange, respectively. We shall later also 
refer to the coupling constants a'^ ^ 0.24 at the hard scale mi,, 
and a" - 0.48 at the ultra-soft scale m^a^. 

Although the string contribution to the potential vanishes for 
light quarkonia in the QGP above it has to be considered at 
T > Tc for heavy quarkonia that remain initially confined and 
are therefore not in thermal equilibrium with the plasma. Hence 
we maintain the string contribution in an approximate solution 
of the gluodissociation problem ifTSll . 




Figure 1: (color online) Radial wave functions of the T(15'), (25), (3S) states 
(solid, dotted, dashed curves, respectively) calculated in the complex screened 
potential eq.(l) for temperatures T = MeV (bottom) and 170 MeV (top) 
with effective coupling constant o-f// ^ (4/3)al = 0.49, and string tension 
cr = 0.192 GeV-. The rms radii < r >''^ of the 2S and, in particular, 35 state 
strongly dependend on temperature T. whereas the ground state remains nearly 
unchanged. 



The string tension of heavy quarkonia decreases with in- 
creasing temperature T in the quark-gluon medium. The 
screened potential including the imaginary part that accounts 
for colhsional damping can be written as ITtI flSl [191 l20l 



V(r, T) = crro [l 
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with roiT) the Debye radius, r^' 



T[4-7Ta'l(2N, + Nf)/6Y'^. 
The number of colors is Nc = 3, the number of flavors in the 
QGP taken as Nf = 3. Because of the inverse proportionality 
of the minimum screening radius that permits a bound state to 
the heavy-quark mass |[T1, it is much more difficult to dissolve 
the T( 1 5 ) in the quark-gluon plasma through screening than the 
JI^{\S). We have calculated the wave functions of the individ- 
ual states, as shown in Fig.[T]for the 15 - 35 states at T = and 
200 MeV. They are almost independent on temperature for the 
ground state. For the 2S state, there is an increase of the rms 
radius from 0.50 to 0.72 fm, whereas for the 35 state the rms 
radius increases from 0.73 to 1.61 fm. 



Due to the high temperature and ensuing large thermal gluon 
density reached at LHC energies in the midrapidity region, the 
most important processes next to screening that lead to a sup- 
pression of Upsilons at LHC are gluodissociation [13 1, and col- 
hsional damping through an imaginary-valued contribution to 
the potential 1 20l [T9l 121112211 (sometimes referred to as Landau 
damping of the exchanged gluon |23 |). 

Whereas gluodissociation of a specific state occurs in a lim- 
ited region of the gluon energy due to the singlet-to-octet ma- 
trix element involved, collisional damping rises with increas- 
ing gluon energy. Hence we calculate the gluodissociation and 
damping cross sections for the y{nS ),Xhi)-P), and^i(2P) states 
as functions of the initial impact-parameter dependent tempera- 
ture in the quark-gluon plasma, and evaluate the time evolution 
in longitudinal and transverse space with collective velocities 
of 0.9c and 0.6c, respectively. Collisional damping is com- 
puted from a numerical solution of the Schrodinger equation 
for a complex potential (see also['24l), with an imaginary part 
as in |20 |. Details of the calculation will be given in 1,25]. 

The leading-order dissociation cross section of the QQ states 
through £1 absorption of a single gluon had been derived by 
Bhanot and Peskin (BP) 1261. Modifying the BP approach 
|[T3l to approximately account for the confining string contribu- 
tion, we use the singlet wave functions computed with Eq.(l). 
Inserting a complete set of eigenstates of the adjoint (octet) 
Hamiltonian -Ajmi, + a"/(6r) with eigenvalues k^/mi, to cal- 
culate the dissociation cross sections of the T(15,25,35) and 
the Xb{^P,2P) states [27], we obtain 
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for the singlet radial wave functions ^^q('") of the b quark, and 
the adjoint octet wave functions §°j('"). The binding energy of 
the nS state is e„, and the 5 function accounts for energy con- 
servation, Irrib - E - e„. 

For vanishing string tension cr — > and the corresponding 
values of the binding energy e„, a pure Coulomb 15 wave func- 
tion, and a simplification in the octet wave function, this ex- 
pression reduces to the result in f2^. 

We obtain new results for the 25 and 35 states from 
eqs. (2), (3). We also calculate the cross sections for th&Xbi^P) 
Xh{2P) states. The gluodissociation cross sections result- 
ing from eqs. (|2]l,(3) including the effect of screening for finite 
string tension are shown in Fig. [2] for the 15 and 25 states. 

One should be prepared to expect modifications in the cross 
section values of the five states from next-to-leading order 
(NLO) contributions ll28l . where a gluon appears in the final 
state in addition to the b and b quarks, and hence, the phase 
space is larger than in leading order (LO). However, in |29| it 
was shown that the quasi-free process that corresponds to NLO 
is less important than LO for temperatures T >270 MeV. 
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Figure 2: (color online) Gluodissociation cross sections crjj^,(nS ) in mb (Ihs 
scale) of the T(15) and T(25) states calculated using the screened wave func- 
tions calculated from the complex potential eq. (1) for temperatures T = 170 
(solid curves) and 250 MeV (dotted curves) as functions of the gluon energy 
Eg. The thermal gluon distribution (rhs scale, solid curve for T = 170 MeV, 
dotted for 250 MeV) is used to obtain the thermally averaged gluodissociation 
cross sections. 



Whereas the heavy quarkonium is not in thermal equilibrium 
with the QGP, it is reasonable to assume that the medium it- 
self is thermalized due to the short equilibration time of about 
0.6 fm/c 1 8 1, at least in the transverse direction. Hence, we in- 
tegrate the gluodissociation cross sections for the T(nS) and 
Xb(nP) states over the gluon momenta p, weighted with the 
Bose-Einstein distribution function of gluons at temperature 
T to obtain the average dissociation cross sections for the nS 
states 
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with E{p) = (p^ + m^)'^^, the gluon degeneracy g,/=16, and 
the gluon density as the integral over the distribution function, 
ng = gdT^((3)/7r^ for nig = 0. Values for the thermal gluon 
density at temperatures 170, 200, 300 and 400 MeV and nig - 
are ng - 1.25, 2.03, 6.85 and 16.23 fm"^, respectively. The 
distribution function is shown in Fig.l2](rhs scale). 

The on-shell gluon energy (p^ + ni^^^^ is usually calculated 
assuming vanishing gluon mass nig = 0. The effect of a finite 
effective gluon mass has been investigated in |fT3l . CoUisional 
damping is added as discussed in |25 1. 

The dissociation widths r{nS,nP) of the nS,nP states are 
then obtained by multiplying the average dissociation cross sec- 
tions with the gluon density. The total widths are obtained by 
adding the corresponding damping widths. In a dynamical cal- 
culation of the fireball evolution ||251 with a functional form 
of the temperature dependence proportional to volume '^'^ (ac- 
cording to the evolution of the energy density), we obtain pre- 
liminary suppression factors R{nS,b) (and analogously for the 
XbinP) states) prior to feed-down at impact parameter b. At 
each impact parameter, the dissociation stops once - 170 
MeV is reached. Related calculations have been performed in 



We then consider the subsequent radiative and hadronic feed- 
down cascade from the T(25,35) and Xh states to the T(15) 
state. In particular, we consider populations estimated from 
the CDF feed-down results |15|. With decay rates for the nS 
states from the particle data group, we calculate a decay cascade 
that matches the final populations measured by CMS for pp at 
2.76 TeV UJJ, and thus provides initial populations which we 
use for the PbPb in-medium calculation at the same energy. 

Following the consideration of screening, gluodissociation 
and damping of the five states, we calculate the radiative feed- 
down cascade in the medium for those states which have sur- 
vived the strong-interaction processes at a given impact param- 
eter b, to obtain the final yields in the presence of the QGP. 

For an T(15) formation time of tf -0.1 fm/c and a quark- 
gluon plasma lifetime of Iqcp - 8 fm/c, our result RAAi^S) ^ 
0.44 is still consistent with the experimental value observed 
by CMS in Q [lO), Raa(IS) = 0.62 ± 0.11(stat)+0.10(sys) 
in minimum-bias PbPb collisions, see also Fig. 3 (Ihs). With 
these values for Y( 15 ) formation time and qgp lifetime, the cen- 
tral temperature at formation time is 670 MeV. The preliminary 
CMS result from the 201 1 run for the ground state suppression 
is 0.53 ||12„|). The suppression factor of the ground state may, 
however, be further reduced by cold nuclear matter effects such 
as gluon shadowing and nuclear absorption which we have not 
considered here. 

For the population ratio of the excited 2S state relative to 
the IS ground state in minimum-bias 2.76 TeV PbPb collsions 
we obtain for an T(15') formation time of tf - Q.l fm/c and a 
quark-gluon plasma lifetime of tgcp - 8 fm/c a model result of 
r(2S)/r(lS)pbPb = 0.35 + 0.15/ - 0.13. Here the theoretical 
error bars arise from uncertainties in the experimental input data 
that we use in the calculation, such as the populations of thexb 
states. The formation times of the excited states - which have 
much larger rms radii - are shorter than the one for the ground 
state. Here we take 0.05 fm/c for the 2S and IP states, and 0.02 
fm/c for the 3S and 2P states. 

CMS has measured a smaller value T(2S)/T(lS)pbpb = 
0.12±0.03(stat)+0.01(sys) (Ha, and an upper hmit of 0.07 for 
the ratio 35/15. Hence, our result (see Fig. 3) leaves room 
for additional suppression mechanisms in PbPb collisions, in 
particular, regarding the excited Y(25) state. 

The expected physical effect, namely, rising dissociation with 
rising temperature, is born out in our approach through the 
combination of screening, gluodissociation, collisional damp- 
ing and feed-down, even though the thermally averaged glu- 
odissociation cross sections first rise and then fall with increas- 
ing temperature for the T{nS) and XbinP) states. 

To conclude, we have calculated the gluodissociation, col- 
lisional damping and screening of T{nS) and Xh(nP) states at 
LHC energies, plus the subsequent radiative feed-down via the 
T(25,35) andxb states. The weakly bound 35 state dissolves 
due to screening already at temperatures T > 180 MeV which 
are close to the critical value. For 25 relative to the 15 state 
we find a substantial suppression due to screening, gluodissoci- 
ation, damping and feed-down that is, however, not as large as 
the value reported by CMS lfT2l . Hence, it allows for additional 
suppression mechanisms of the excited state. 
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Figure 3: (color online) Suppression factor Raa for the T(15) state (Ihs), ratio T(2S)/T(15) (middle), and T(3S)/T(15) (rhs) calculated in the present work for 
2.76 TeV PbPb minimum-bias collisions from screening, gluodissociation, collisional damping and feed-down for a 15 formation time tf = QA fm/c as functions 
of the quark-gluon plasma lifetimes tQcp. The corresponding CMS minimum bias results are shown from the 2010 run 1 11 1 for Raa(^S) (Ihs, solid line with dotted 
statistical and systematic uncertainties), and preliminary results from the 201 1 run 1 12 1 (Ihs, dashed line). The middle graph for 25 /IS shows our results (bars) above 
the preliminary CMS data from 1 12| (solid line, with experimental uncertainties dotted), leaving room for additional suppression mechanisms, or modifications of 
the formation time scales. The 35/15 calculations (rhs) are compared with the preliminary CMS upper limit (dashed line 1 12|). The estimated theoretical error bars 
account for the uncertainties in the input data that enter our calculations. 



We obtain results for the suppression of the excited T states 
relative to the ground state in PbPb coUsions at LHC energies 
with T(15) formation time 1 fm/c and QGP lifetimes of up to 
1 1 fm/c. Typical central QGP temperatures have a strong time 
dependence. For central collisions, the initial temperature at a 
r(lS) formation time of 0. 1 fm/c is 670 MeV. 

Although screening of the strongly bound 15 ground state is 
negligible, we find that its gluodissociation is sizeable due to 
the strong overlap of the IS gluodissociation cross section with 
the thermal gluon distribution, and also collisional damping is 
relevant. The observed suppression factor Raa{^S) - 0.62 in 
minimum-bias PbPb collisions ifTTI (0.53 in the preliminary 
2011 data [12|) is mainly due to both direct gluodissociation 
and damping of the 15 state, and to the melting and gluodisso- 
ciation of the Xbi^P) ^ndxb(2.P) states which partially feed the 
15 state in pp, pp and e^e^ collisions. 

The preliminary CMS data from the 2011 run that are 
available now [12] provide also the centrality dependence of 
the suppression factors for both 15 and 25 states. We shall 
soon compare with these data 1,25 J . 

Acknowledgments 

This work has been supported by the ExtreMe Matter Insti- 
tute EMMI, and IMPRS-PTFS Heidelberg. 



[10] S. Chatrchyan, et al., CMS-PAS-HIN-10-006. 

[1 1] S. Chatrchyan, et al., Phys.Rev.Lett. 107 (201 1) 052302. 

[12] T. Dahms, et al., Proc. Hard Probes (2012) in press. 

[13] F. Brezinski, G. Wolschin, Phys. Lett. B 707 (2012) 534-538. 

[14] Y. Mehtar-Tani, G. Wolschin, Phys. Rev. Lett. 102 (2009) 182301. 

[15] T. Affolder, et al., Phys.Rev.Lett. 84 (2000) 2094. 

[16] E. Eichten, K. Gottfried, T. Kinoshita, J. Kogut, K. D. Lane, T.-M. Yan, 

Phys.Rev.Lett. 34 (1975) 369. 

[17] S. Jacobs, M. G. Olsson, C. Suchyta, Phys. Rev. D 33 (1986) 3338. 

[18] F Karsch, M. Mehr, H. Satz, Z.Phys.C 37 (1988) 617. 

[19] A. Beraudo, J. Blaizot, C. Ratti, Nucl. Phys. A 806 (2008) 312. 

[20] M. Laine, O. Philipsen, P Romatschke, M. Tassler, JHEP 03 (2007) 054. 

[21] N. Brambilla, J. Ghiglieri, A. Vairo, Phys. Rev. D 78 (2008) 014017. 

[22] N. Brambilla, et al., Eur.Phys.J.C 71 (2011) 1534. 

[23] M. Strickland, Phys. Rev. Lett. 107 (2011) 132301. 

[24] J. Wahner, Heidelberg BSc thesis (2012). 

[25] F. Nendzig, G. Wolschin, to be submitted to PRC. 

[26] G. Bhanot, M. E. Peskin, Nuclear Physics B 156 (1979) 391. 

[27] F. Nendzig, Heidelberg PhD thesis (in preparation). 

[28] T. Song, S. H. Lee, Phys. Rev. C 72 (2005) 034002. 

[29] L. Grandchamp, R. Rapp, Phys. Lett. B 523 (2001) 60. 

[30] T. Song, K. C. Han, C. M. Ko, Phys. Rev. C 85 (2012) 014902. 



References 

[1] T. Matsui, H. Satz, Phys. Lett. B178 (1986) 416. 

[2] L. Kluberg, H. Satz, arXiv:0901.3831. 

[3] D. Kharzeev, J. Phys. GiNucl. Part. Phys. 34 (2007) S445. 

[4] E. Atomssa, et al., Nucl.Phys.A. 830 (2009) 331c. 

[5] B. K. Patra, D. K. Srivastava, Phys. Lett. B 505 (2001) 1 13. 

[6] G. Martinez Garcia, et al., J. Phys. G: Nucl. Part. Phys. 38 (201 1) 124034. 

[7] C. Silvestre, et al., J. Phys. G: Nucl. Part. Phys. 38 (201 1) 124033. 

[8] C.-Y. Wong, Phys.Rev.C 72 (2005) 034906. 

[9] H. Masui, et al., J. Phys. G: Nucl. Part. Phys. 38 (201 1) 124002. 



4 



